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by analysis of experimental and numerical data, has its origin in the binding properties of the ionic 
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I. INTRODUCTION 

In a variety of applications in physics it is important to 
understand how is the normal Hydrogen spectrum mod- 
ified if at small distances the Coulomb law is replaced by 
a central short-ranged potential. An important aspect of 
this problem is the existence of two length scales - the 
Bohr radius of the Coulomb field as and the range of 
action of short-range forces ro . 

For example, in hadronic atoms formed by charged par- 
ticles and antiparticles the large distance Coulomb at- 
traction gives its way at short distances to nuclear forces 
whose range ro is significantly smaller than as 0]. 

In condensed matter physics a similar problem is that 
of the energy spectrum of the Wannier-Mott exciton Q. 
When in a semiconductor an electron is excited into the 
conduction band, a bound state with a hole left in the va- 
lence band can form. Due to the large dielectric constant 
of the medium the electron and the hole in the exciton 
are spatially well-separated. Therefore the electron-hole 
interaction is a Coulomb attraction modified at short dis- 
tances. In this context as can exceed many times ro 
which is of the order of the Hydrogen Bohr radius. 

Zel'dovich was apparently the first to recognize that in 
the limit ro -C as the spectrum of the distorted Coulomb 
problem is peculiar Q. Since the centrifugal barrier de- 
creases the probability of particle penetration in the re- 
gion of small distances r, the effect of the short-range 
potential is strongest for the states of zero angular mo- 
mentum. In this case the radial motion of a particle of 
mass m and energy E in a central potential U(r) is de- 
scribed by the one-dimensional Schrodinger equation [J] 



fx 

dr 2 



+ — (E-U(r)) X = 



(1) 



where xi r )/ r ls the radial wave function. Zel'dovich 
chose U(r) — —h 2 /maBT for r > ro and U(r) — U s (r) 



otherwise and demonstrated that as long as the short- 
range potential U s (r) is not resonant, its effect is weak. 
If, on the other hand, U s (r) has a low-energy scatter- 
ing resonance, a drastic reconstruction of the spectrum 
takes place. Using the example of the square well of 
depth Uq he stated that as the dimensionless coupling 
constant w ~ mr^Us/h 2 increases, the spectrum of the 
problem E n {w) evolves in a fashion resembling a sharp 
decreasing staircase. The steps are located at critical val- 
ues of w at which bound states occur in U s (r) only. As 
w goes through the first threshold, the Coulomb levels 
E n (n > 2) quickly fall to E n -\ while the ground state 
E\ rapidly drops downward. The relative width of the 
region where the spectrum reconstruction takes place, 
Aw/w ~ ro/as -C 1, is narrow, and qualitatively the 
same pattern repeats itself upon passing through every 
subsequent resonance. 

A similar spectral behavior has been found by Popov 
Q in his analysis of the Dirac equation for an electron in 
a field of the bare nucleus of charge Ze with Z > 137. 

The Zel'dovich effect has been re-discovered in the 
spectra of hadronic atoms, and its generality has been 
demonstrated for any interaction with two widely differ- 
ent spatial scales 0. 

Various aspects of the spectrum reconstruction have 
been investigated by Popov and collaborators Q. Their 
study was motivated by then existing experimental evi- 
dence of the large ls-level shift in the proton-antipropton 
atom which was naturally linked to the Zel'dovich effect. 
Later it became clear that the experimental level shifts 
are small and the interest in the phenomenon declined. 

As far as we know, at this time there is no documented 
experimental evidence of the Zel'dovich effect. This is 
not surprising because the spectrum reconstruction takes 
place in a narrow range of parameters in the vicinity of 
low-energy resonances. However a given experimental 
system is unlikely to be near resonance. A systematic 
search for the Zel'dovich effect would consist in looking 
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for spectral changes in response to tuning of the cen- 
tral part of the potential which is often impossible - the 
strength of the nuclear force cannot be changed in the 
laboratory. 

Recently Karnakov and Popov y| pointed out that the 
Zel'dovich spectrum reconstruction takes place for a Hy- 
drogen atom as a function of the external magnetic field 
thus providing an example of a system where a systematic 
search for the effect might be possible. Although the phe- 
nomenon is observable in numerical studies, direct exper- 
imental evidence is lacking and may only come from as- 
trophysical observations as the pertinent magnetic fields 
are comparable to those on the surface of a neutron star. 

The goal of this paper is to demonstrate that evolu- 
tion of the Rydberg spectra of ordinary atoms along the 
Periodic Table provides direct evidence of the Zel'dovich 
effect. Since the condition ro <C as does not hold in 
atomic systems, the way the phenomenon manifests it- 
self is less dramatic - we will show that it can be seen as 
a systematic periodic spectral modulation as a function 
of the cubic root of atomic number Z . 

It is known that for a highly excited s electron of a Ry- 
dberg atom the effect of polarization of the ionic core is 
negligible compared to that of the wave function penetra- 
tion in the central region of the atom [flj . Therefore the 
electron dynamics can be adequately described by Eq.flJ 
where the effective central field U (r) at large distances is 
that of a positively charged ion of charge e. On the other 
hand, starting from distances of the order of the size of 
the ionic core ro ~ as the field felt by the electron begins 
to deviate from the —e 2 /r form on-average decreasing, as 
r —>■ 0, to —Ze /r. By increasing Z along the Periodic 
Table Nature systematically deepens the inner part of 
the potential leaving the outer —e 2 /r tail intact. Thus 
by analyzing the Rydberg spectra as a function of atomic 
number Z it may be possible to correlate them with the 
binding properties of the ionic core which will constitute 
evidence of the Zel'dovich effect. 

In atomic physics the motion of an electron in the field 
of a residual atomic ion has been studied in the past. 
Approximating the potential of the ionic core by that of 
the Thomas-Fermi or Thomas-Fermi-Dirac theories Lat- 
ter ^3] computed numerically the single-electron term 
values from Is to 7d for all atoms. His ns spectra as a 
function of atomic number Z for largest n studied clearly 
show modulations on a decreasing energy curve. It is 
well-known that the large n atomic spectra are described 
by the Rydberg formula 

h2 1 

En — ~7i T7 \9 ( 2 ) 

2ma B (n — (i) z 

where fi is the quantum defect which in the limit n — > oo 
does not depend on n. Latter's results imply that the de- 
pendence of the quantum defect ponZ has modulations 
superimposed on an increasing curve. 

The fi(Z) dependence has been numerically computed 
by Manson [ll| and by Fano, Theodosiou and Dehmer 



[l2T | who used the Hartree-Slater model to approx- 
imate the potential of the ionic core of the atom. Al- 
though the periodic variations of fi(Z) are strongly ob- 
scured by the shell effects (included in the Hartree-Slater 
model), Fano, Theodosiou and Dehmer argued that they 
are there and that there is a correlation between the 
location of radial nodes of the function x from Eq.flJ 
near r$ and the slope of the fi{Z) dependence. In view 
of the oscillation theorem |3| the nodal structure of the 
function \ is intimately related to the binding properties 
which suggests that systematic periodic variations of Ry- 
dberg spectra as function of Z might be related to the 
Zel'dovich effect. 

In order to demonstrate that this connection is correct 
below we compute the upper part of the spectrum of the 
modified Coulomb problem not assuming that ro <C ag. 
We show that the staircase reconstruction taking place 
for ro < «b and the spectral modulations for ro ~ as 
are different limiting cases of the same phenomenon - 
sensitivity to the binding properties of the inner part of 
the potential which we continue to call the Zel'dovich 
effect. We also compare our results for fi(Z) with avail- 
able experimental and numerical data to show that the 
phenomenon is observable. 

The organization of this paper is as follows. In Section 
II we provide a short derivation of the Rydberg formula 
<j2J and arrive at the expression for the quantum defect 
in terms of the dimensionless range of the inner potential 
and its scattering length. This general result is further 
analyzed in the ^ C «b limit and the main features 
of the Zel'dovich spectral reconstruction are recovered 
(Section IIA). In Section IIB we establish a relationship 
between the Zel'dovich effect and Levinson's theorem of 
quantum mechanics. This is followed (Section IIC) by 
the analysis of the opposite ro 3> as limit where we 
demonstrate that the Zel'dovich effect manifests itself in 
the form of a spectral modulation whose origin still lies in 
the binding properties of the inner potential U s (r). These 
general findings are illustrated in Section IID where we 
use the exactly-solvable example of the rectangular well 
as a model for the inner potential. In Section HE we 
observe that only a treatment more accurate than semi- 
classical can capture the Zel'dovich effect. 

Section III focuses on the computation of the system- 
atic quantum defect of the Rydberg electron as a function 
of atomic number Z. First (Section IIIA), for the inner 
potential having an attractive Coulombic singularity at 
the origin, we derive a semiclassical expression for the 
quantum defect and show that it is equal to the number 
of de Broglie's half- waves fitting inside the inner potential 
minus a contribution proportional to (ro/as) 1 / 2 . Going 
beyond the semiclassical approximation we also demon- 
strate that the Zel'dovich modulation of the quantum de- 
fect is a periodic function of the number of de Broglie's 
half-waves fitting inside the ionic core of the atom. This 
is followed by an explicit calculation based on Latter's 
model of the ionic core First, the semiclassical quan- 
tum defect is calculated as a function of Z 1 / 3 (Section 
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IIIB). Then (in Section IIIC) a full computation cap- 
turing the Zel'dovich effect is performed. An important 
ingredient here is an approximate calculation of the scat- 
tering length of the ionic core of the atom. Both the 
scattering length and the related Zel'dovich modulation 
of the quantum defect turned out to be nearly periodic 
functions of Z 1 / 3 . 

In Section IV the results of our systematic calcula- 
tion are compared with experimental and numerical data. 
First, we observe that the bulk of the quantum defect val- 
ues is well-captured semiclassically. Then (Section IVA) 
we demonstrate that the gross features of the deviation 
away from semiclassics are due to the effects of the shell 
structure. This is done by establishing and demonstrat- 
ing a correlation between the variation of the radius of 
the ionic core of the atom and corresponding variation 
of the quantum defect. Finally, in Section IVB a Fourier 
analysis of the quantum defect variation with Z x l 3 is con- 
ducted which singles out the Zel'dovich effect. As a by- 
product we also find a Z 1 ! 3 periodic contribution coming 
from the shell effects. 

We conclude (Section V) by outlining our main result 
and directions of future work. 



II. DISTORTED COULOMB PROBLEM AND 
QUANTUM DEFECT 

We will be interested in the low energy bound states 
with the classical turning point being far away from the 
boundary of the central region, i. e. h 2 /ma,B\E\ ^> ro. 
Then the quickest way to derive the spectrum is via semi- 
classical arguments derived from those given by Migdal 

For tq < r < h 2 /maB\E\ the semiclassical solution to 
Eq.Q can be written in two equivalent forms: 



1 
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(3) 



where p = (—2m\E\ + 2h 2 /ras) 1 ^ 2 is the momentum. 
The first representation in Eq.<EJ is the standard result 
with the phase of 7r/4 improving on the deficiency of 
the semiclassical approximation near the classical turning 
point, while the yet undetermined phase a in the second 
representation in Eq.(EJ both corrects for the failure of 
the semiclassical approximation in a Coulomb field at dis- 
tances r <clb and accounts for the short-range potential 
U s (r). 

For ro < r <C h 2 jmaB\E\ the Schrodinger equation 
simplifies to 



fx 

dr 2 



ran 



-x = o 



(4) 



and can be exactly solved: 

X oc r 1 ' 2 (,hW%r/a B ) - Y^Zr / a B )t & n&) (5) 

where J v {x) and Y v {x) are the order v Bessel functions of 
the first and second kind respectively The solution 
© is a linear combination of the regular J\ (0) = and ir- 
regular 1^(0) = oo Coulomb functions of zero energy, and 
for the purely Coulomb problem, U s (r) = —Yc 2 jmasr, 
one has to recover tan 5 = 0. 

For as <C r <C h 2 /ma,B\E\ the semiclassical approx- 
imation is accurate, and the second rep resentat ion of 
Eq.© yields x r 1 / A sin( i/S^ob — \Z&r a /aB + a). 
On the other hand, the r ^> as limit of © is x oc 
r 1 / 4 sm(y/8r /as — tt/4 + 5) which determines a in (01 
to be y/^ r o/ a B — tt/4 + S. It also implies that 5 in © 
is the zero-energy phase shift due to the small-distance 
deviation of the potential from the Coulomb form. 

The energy spectrum can be found from the require- 
ment that the semiclassical expressions (0 coincide. 
Combined with a = y/ 8ro /as — tt/4 + 5 this gives the 
quantization rule 



1 

h 



h 2 /ma B \E\ 



pdr 



x 



where dimensionless parameter 



X 



(6) 



(7) 



measures the range of the short-range forces. Calculating 
the integral we arrive at Eq.(|2j) with S = irfi which is the 
statement of Seaton's theorem 01 relating the quantum 
defect to the zero-energy phase shift. 

The range of applicability of Eq.(|2l, n — fi ^ 
(ro/as) 1 / 2 ~ x , follows from the condition \E\ <C 
h 2 /maBro which also implies that in order to calculate 
the quantum defect entering the spectrum i ffi l. we only 
need to match (0 with its zero energy counterpart at 
r < r Q . 

We proceed by computing h = [dlnx/dlnr] r ^ ro+ o ; 
the logarithmic derivative of the function (JHJ evaluated 
at the boundary of the inner region: 



h = 



x Jq(x ) - Yo(x )tan7r^ 
2 Ji(iro) ~ Yi(xq) tan7r/x 



(8) 



where we used S = 7r/i. The quantum defect [i is deter- 
mined by equating (0 to h s = [dlnx/dln v\ r — > ro — o which 
can be found by solving the E = Schrodinger equation 
for r < ?' with U(r) — U s (r): 



a-X 
dr 2 



2m 
H 2 



U s (r) X = 



(9) 



The parameter h s can be equivalently expressed in terms 
of the scattering length corresponding to the inner po- 
tential only. Indeed for motion in a short-range potential 
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27ta = 




Figure 1: Quantum defect for U s (r) = as a function of the 
range parameter xo, Eq.JJJl, and its io > 1 limit, fi(0, xq) = 
3/4 — xq/tt (shown in gray scale). 



Figure 2: Graphical solution of En. 11211 : xo = 1/30 has been 
used to construct the graph. The quantum defect [i is given 
by the intersections of the right-hand side of l|12fl with the 
line of constant as/27ra s . 



the scattering length a s is defined from the asymptotic 
r — > oo behavior \(r) oc 1 — r/a s of the solution to yj). 
For a potential well identically vanishing for r > r , this 
is also the exact behavior outside the well with the im- 
plication that |3J, [3| 



_ l' dlnx(r -> r - 0) 
rflnr 



'•() 



(10) 



Then substituting h = h s in Eq.® and using (ff]~0f) we 
arrive at the formula for the quantum defect 



2x 1 Ji(x ) + (a a /rQ - I)Jq(xo) , . 

tan7r/i = (11 

1 Y 1 (x ) + (a s /r - l)Y (x Q ) 



If the short-distance potential is selected in the form 
U s (r) — —ti 2 /masr, i. e. we have the ordinary Coulomb 
problem in the whole space, the quantum defect /i enter- 
ing the Rydberg formula Q must vanish identically. It is 
straightforward to verify that this is indeed the c ase: th e 
E = inner r < tq solution to QJ, x oc r 1 / 2 Ji(y/&r /<xb), 
leads to the expression for the scattering length nullify- 
ing the numerator of l(TT|l . This argument defines the 
zero of the quantum defect and implies that [i is neces- 
sarily positive if for all r < rg the inner potential U s (r) is 
more attractive than the Coulomb potential —H 2 /maBr; 
otherwise the quantum defect is negative. For example, 
for U s (r) — the quantum defect /i(0,xo) is a neg- 
ative monotonically decreasing function of xq such as 
ji = — Xq/32 for io 1, and /i = 3/4 — x /ir in the 
opposite x 3> 1 limit. The fx(Q,x ) dependence as well 
as its Xq 3> 1 limit are shown in Fig. 1. 



A. Zel'dovich effect in the ro as limit 

For xo <C 1 Ea. llll|l simplifies to a form accumulating 
the physics of the Zel'dovich effect: 



as 
2nae 



nx 2 M 



2 2 
cot7r/i In (12) 

7T 7X 



where In 7 — 0.5772 is Euler's constant. Terms of higher 
order in Xq which for U„(r) — lead to small negative 
values of the quantum defect are neglected in J12II . 

We verified that Eo. l(T2)l matches the upper portion of 
the ns spectrum which for ro <C as is known in closed 
form for any n Q. We also note that with some effort 
Eq. i112ll can be deduced from the expression for the phase 
shift of the proton-proton scattering given by Landau and 
Smorodinskii 0|: in their formula we have to (i) reverse 
the sign of the Bohr radius, (ii) take the limit of zero 
energy, and (iii) employ Seaton's theorem Q S = nfi. 

Fig. 2 shows the dependence of aB/2na s on [i given by 
Eci, l|12|l : its inverse /i(as/27ra s ) is a multivalued function 
consisting of a series of increasing step-like curves sand- 
wiched between nearest non-negative integers. The slope 
of fi(aB/2na s ) is small everywhere except for the vicinity 
of half-integer 

Since typically the central potential U s (r) is not res- 
onant, \h s \ in ltl2|l is not small. Then the magnitude 
of the scattering length is of the order of the size of 
the inner well, \a 8 \ ~ ro and as/27r|a s | ~ l/^o ls sig- 
nificantly larger than the last term in I112L This im- 
plies that the quantum defect is very close to an integer, 
/j = — 2a s /a,B (mod 1), with \a s \/as — ro/sB <§C 1. This 
conclusion is in quantitative agreement with the results 
of perturbation theory in a s /as when the deviation from 
the Bohr Hydrogen formula is small 0, H|- It is ap- 
plicable to an attractive non-resonant well of arbitrary 
strength; for weak U s (r) which cannot support a bound 
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state we have /i = — 2a s /aB > represented by the 
leftmost intersection in Fig. 2. We also note that the 
spectrum is exactly Hydrogenic if the scattering length 
is zero which can be viewed as an analog of the Ramsauer 
effect [19]: in the present context it refers to a resonant 
phenomenon when the distortion of the Coulomb poten- 
tial at small distances is invisible to the low-energy bound 
(or incident) particle. 

Exactly at half-integer fi the scattering length is neg- 
ative with the magnitude \a s \ — (as/4) In -1 (2/7x0) = 
2tqXq 2 In -1 (2/7x0) significantly exceeding the size of the 
central region tq. This implies that the slope of the 
/j,(aB/2Tra s ) dependence is largest when U s (r) itself is 
almost resonant so that it supports a low-energy virtual 
state. At the point of the steepest slope we also have 
[dx/dr) r=ro cx h s ~ (xq/2) ln(2/7Xo) <C 1. Since this is 
practically zero, one can equivalently say that the slope of 
the n{dB I^kcls) dependence is largest when the antinode 
of the function x m Eq.Q occurs at the boundary of the 
inner region r$. This criterion resembles that given by 
Fano, Theodosiou and Dehmer 0] for the dependence 
of the quantum defect [i on atomic number Z. We note 
however, that for a Rydberg atom the size of the residual 
ion does not satisfy the condition r <C as; this issue is 
further addressed below. 

If for all r the central well is attractive, its effect can 
be quantified by a single dimensionless coupling constant 
w ~ mrl\U s \/h 2 > where \U S \ has a meaning of the 
characteristic depth of the well. Then the inverse scatter- 
ing length a J 1 is known to be a monotonically increasing 
function of w Q - an a s dependence shown in Fig. 3 in 
gray scale is typical and may help illustrate the argument 
given below. 

The step-like features of the function fi(aB/2na s ) are 
amplified in the fi(w) dependence. Indeed, for w<l the 
scattering length a s is very small and negative. Then the 
line of constant aB/2ira s in Fig. 2 lies at very large nega- 
tive values, the quantum defect satisfies [i = —2a s /aB <C 
1 , and the deviation from the normal Hydrogen spectrum 
is small. As the well deepens, the coupling constant w 
increases, the scattering length becomes more negative 
and the horizontal line of constant as/27ra s moves up- 
ward. However as long as the well remains non-resonant, 
the quantum defect fi will only grow very little. The 
strongest increase of /Lt(u>) in response to deepening of the 
well (and thus the largest deviation from the Bohr Hydro- 
gen formula) occurs when the scattering length reaches a 
very large negative value a s = —2toXq 2 In -1 (2/jxq). For 
10 C 1 this takes place very close to a threshold value 
of the coupling constant w when the first bound state is 
about to appear in U s (r). The relative width of the re- 
construction region Aw/w centered around n — 1/2 thus 
can be estimated from the scaling behavior of the scat- 
tering length near the threshold a s ~ r w/Aw and the 
condition or ~ \a s \. This leads to the original result of 
Zel'dovich Q Aw/w ~r /aB- 

As the coupling constant w increases through the first 
binding threshold, the inverse scattering length changes 
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Figure 3: Evolution of the Zel'dovich effect for the rectangular 
well of radius ro and depth Uo for a series of range parameters 
xo, Eq.Q, manifested in the dependences of the quantum 



defect (ion(= (8mUor Q /n h 
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The correlation 



with the binding properties of the well is seen from the plot of 
the reduced scattering length a a /ro (gray scale). Shown are 
also the values of the zero-energy phase shift S = Tifi relating 
the Zel'dovich effect to Levinson's theorem. 



sign, and the line of constant as/27ra s in Fig. 2 enters 
the region of positive values. After passing through the 
reconstruction region, the positive scattering length de- 
creases in magnitude, for a s /as <C 1 the quantum defect 
is close to unity, /1 = 1 — 2a s /as, and the deviation from 
the normal Hydrogen spectrum is again small. In the re- 
gion a s ~ ro the scattering length does not vary strongly 
with the depth of the well, and one can say that the 
slope of the /j,(w) dependence will be minimal when the 
node of the function x m Eq.Q is near the boundary 
of the central region r which parallels the criterion of 
Fano, Theodosiou and Dehmer Upon further in- 

crease of the coupling constant w, the scattering length 
gets smaller and the line of constant aB/2ira s in Fig. 1 
enters the region of very large positive values becoming 
infinite at a s = 0. 

To summarize, as a s (w) goes through one complete 
cycle decreasing from zero, passing through the binding 
resonance, and then approaching zero from above, the 
quantum defect fi(w) increases from zero to unity in a 
staircase fashion: it is mostly zero or unity except for 
the narrow region Aw/w ~ r^/as <C 1, [i — 1/2 near 
the first binding threshold of U s (r). Combined with the 
Rydberg formula 1J2J this implies that the Coulomb levels 
E n quickly fall to -E„_i which constitutes the essence of 
the Zel'dovich effect Q. 

As the coupling constant w continues to increase away 
from a s (w) = 0, the next cycle, 1 < fi(w) < 2, begins and 
qualitatively same pattern repeats itself. This remains 
true for every subsequent cycle with fi(w) sandwiched 
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between nearest integers. Overall the quantum defect is 
an increasing function of w having the form of a staircase 
with practically integer plateaus and sharp steps located 
at half-integer fi. The steps correspond to the presence 
of the low-energy scattering resonances in U s (r) . 

To illustrate this behavior we choose the inner poten- 
tial in the form of a rectangular well of depth Uq whose 
scattering length is given by a s /r = 1 — 2tan(7r£/2)/7r£ 
with dimensionless parameter £ = (SmUorl/iT 2 ^, 2 ) 1 / 2 ~ 
W 1 / 2 quantifying the depth of the well. The scattering 
length diverges at odd values of £ which correspond to 
consecutive occurrences of bound states in the well; the 
respective dependence of a s /r on £ is shown in Fig. 3 
in gray scale. We also plot the dependences of the quan- 
tum defect /i on £ found from the general expression 11 111 
for a series of representative x$. The analysis based on 
Eo. l(T2)l is illustrated by the xq = 1/10 and xo = 1/2 
curves; the latter corresponds to the case of the proton- 
antiproton atom 0|. These dependences have the form of 
staircases with nearly integer plateaus; the steepness of 
the steps where the quantum defect varies by unity and 
the flatness of the plateaus increase as xq gets smaller. 
An inspection reveals that the points of maximal slope 
of somewhat precede the scattering resonances in 
accordance with the analysis given above. This is seen 
most clearly for the £ ~ 1 step of the xq — 1/2 curve. 
Fig. 2 of Zel'dovich's work 3j| has this feature as well. 
From a practical standpoint the steps can be considered 
to coincide with the binding resonances of the well. 

The relative width of the reconstruction region A£/£ 
can be estimated as ro/as£ 2 - Since the threshold values 
£ grow linearly with the number of bound states, then 
for fixed xq the steepness of the steps increases with £ as 
can be seen in Fig. 3. This is merely the consequence of 
the sharpening of the binding resonances. Similarly the 
flatness of the plateaus improves as £ increases, and the 
points of least slope of the /i(£) dependence asymptoti- 
cally approach even values of £. This is where the node 
of the function x m Eq.QJ coincides with the boundary 
of the central region, a s = ro. 

Finally we note that the quantum defect takes on ex- 
actly integer values whenever the scattering length van- 
ishes. 



B. Connection to Levinson's theorem 

There is a deep parallel between the Zel'dovich recon- 
struction of the upper E — > part of the Coulomb spec- 
trum in the r$ <C clb limit and the low-energy scattering 
by a short-range potential well. For a particle of energy 
E = ti 2 k 2 /2m whose wave vector k is small in magni- 
tude, kro <§C 1, scattered by the short-range potential 
U s (r) vanishing for r > tq the scattering length a s can 
be defined 20] through the k — > limit of the relationship 

\/ka s = -cot 8 s (k) (13) 



where S s (k) is the phase shift. Employing Seaton's theo- 
rem S = 7T/i relating the quantum defect to the zero- 
energy phase shift it is straightforward to realize that 
Eas. <ll2H and 111- -St are direct analogs. The Coulomb field 
is characterized by its own length scale, the Bohr ra- 
dius as- Its free particle counterpart entering Eo. ifT3|l 
is the de Broglie wavelength 27r/fc. The range of ap- 
plicability of Eci. lll2H ro < hb parallels the low-energy 
condition kro "C 1 necessary for Ea, l|13Jl to hold. The 
analysis which led to the explanation of the Zel'dovich 
effect can be repeated for Eo. lfl3ll with the conclusion 
that the phase shift S s (k) as a function of the dimen- 
sionless depth w of the scattering well has the form of a 
sharp increasing staircase whose plateaus practically co- 
incide with 5 s {k) = (mod n). The steps where the phase 
shift changes by ir are very narrow, Aw/w ~ kro <C 1, 
and the points of steepest slope are located at 5 s (k) ~ 
7r/2 (mod 7r). In the limit k — > the staircase becomes 
perfect. This can be recognized as Levinson's theorem 
|20l | relating the number of bound states in a well with 
the zero-energy scattering phase shift. We conclude that 
for r <C as the Zel'dovich effect expressed in terms of 
the zero-energy phase shift S is the Coulombic cousin of 
Levinson's theorem jjjjjj. A special case of this correspon- 
dence, the Ramsauer-like recovery of the normal Hydro- 
gen spectrum for a s = 0, w as already mentioned earlier. 
In the limit xq = \J 8ro /as — * Zel'dovich's staircase be- 
comes perfect and identical to Levinson's staircase. This 
can be understood as a result of taking the neutral limit, 
cib — » oo, when the Coulomb part of the binding poten- 
tial U(r) in vanishes. From this viewpoint, Levinson's 
theorem is a consequence of the Zel'dovich effect. To em- 
phasize the connection to Levinson's theorem, in Fig. 3 
we additionally show the zero-energy phase shift 8 = nfi. 

Fig. 3 also demonstrates that as xq increases, the 
staircase //(£) dependence with well-defined steps and 
plateaus evolves into an increasing function with mod- 
ulations: the "plateaus" develop noticeable slope and the 
"steps" acquire a width. Moreover for sufficiently large 
x > 1 the staircase-like appearance seems to emerge 
only for a sufficiently deep well, i. e. large £. Another 
feature is the presence of a negative offset which is a 
growing function of xo . This is due to the fact that for 
U s (r) — the quantum defect is a monotonically decreas- 
ing negative function of xq as shown in Fig. 1. 



C. Zel'dovich effect in the ro 2S> as limit 

In the xq ^> 1 limit Ea Jllll simplifies to the form 

a s 2 . 37T, 

— = — cot(7r^ + .T r) + ! (14) 

r x 4 

allowing model- independent treatment. The analysis of 
Eq. ifl^Ijl is convenient to conduct in terms of the reduced 
quantum defect M — ^ + x /tt — 3/4 whose zero gives the 
xq ^> 1 asymptotic of /i for U s (r) = 0. Ea. ill4|) can be 



M 



Figure 4: Graphical solution of En. ll4|l : xo = 4 has been 
used to construct the graph. The reduced quantum defect 
M = H + xo /ir — 3/4 is given by the intersections of the right- 
hand side of 114t with the line of constant a s /ro. 



Figure 5: Plot of the surface of the relative quantum defect 
A/j,(£,xo) = /i(C, x o) — ^(0,xo) for a rectangular potential 
well with dimensionless range and depth parameters xo and 
C, respectively, according to Ea. lllll . 



investigated in a manner analogous to that of Eq. lf3~2f) ; a 
brief summary is given below. 

Fig. 4 shows the dependence of a s /r on M = fj, + 
x /tt — 3/4 given by En. ltT^jl : its inverse M(a s /r ) is a 
multivalued function consisting of a series of decreasing 
step-like segments sandwiched between nearest integers. 

The magnitude of the slope of M(a s /ro) is smallest at 
integer M which occurs at binding resonances, a 8 = ±oo, 
i. e. when the antinode of the function \ m Eq.lQJ 
coincides with the boundary of the inner region tq. In the 
vicinity of integer M we find M = 2ro/irxoa s (mod 1). 
This translates into an explicit result for the quantum 
defect [i = 3/4 — xq/tt + 2ro/^Xoa s (mod 1) valid in the 
limit xq ^S> 1 and r /xoa s <C 1, thus roughly covering the 
range of \a s | from tq to infinity. In the vicinity of the first 
binding resonance we have [i = 3/4 — xq/tt + 2ro/irxoa s 
which is represented by the leftmost intersection in Fig. 4. 

The magnitude of the slope of the M(a s /ro) depen- 
dence is largest at half-integer M which occurs at a s = ro, 
i. e. when the node of the function \ m Eq.Q coin- 
cides with the boundary of the inner region ro. Since 
the reduced quantum defect M is a decreasing function 
of a s /ro, and the scattering length a s is a decreasing 
function of the well depth w @], then for fixed xq the 
parameter M (and thus the original quantum defect /x) 
is an increasing function of w. 

In contrast to the io « 1 regime, here the step- 
plateau features of the function M(a s /ro) are generally 
suppressed in the M(w) dependence. This is because the 
dependence of the scattering length a s on the depth of 
the well w is weakest in the region a s — ro where the 
M(a s /ro) dependence shows a "step". By the same to- 
ken the "plateaus" acquire a noticeable slope since the 
a s (w) dependence is strongest near the binding reso- 
nance, a s = ±oo, i. e. where the M(a s /r ) dependence 
is weakest. As a result the dependence of the quantum 



defect fi on the depth of the well w is more appropriately 
viewed as consisting of modulations superimposed on an 
increasing curve. These modulations still have their ori- 
gin in the binding properties of the inner potential U s (r). 



D. Rectangular well example 

The analysis of Sections IIA and IIC is illustrated in 
Fig. 5 where using the example of the rectangular well 
and Eq. ifTTTl we plot the surface of the relative quantum 
defect A[i((, xq) = [i((, x ) — /x(0, xo). The point of sub- 
tracting ^t(0,xo) from /i(C,a;o) is to isolate the physics of 
binding from the background /Lt(0,xo) which is a mono- 
tonically decreasing function of xo shown in Fig. 1. The 
peculiar shape of the resulting surface can then be un- 
derstood as follows: 

For x <C 1 the background contribution /i(0, xo) 
is negligible (see Fig. 1), the relative quantum defect 
A/x(C, xq) reduces to /x(£, Xq) which, according to our ear- 
lier analysis, is a staircase function of £ with steps located 
at odd C, i. e. when the bound states occur in the well. 

For xo ^> 1 the background contribution fi(0,xo) is 
3/4 — xq/tt, and we find that Afi(£,xo) — /i(£,£o) — 
n(0,xo) = fi((,x ) + x /n — 3/4 = M which, accord- 
ing to Fig. 4, is a decreasing staircase function of a s /ro- 
For not very deep well the dependence on ( has a form 
of a rounded staircase with "steps" and "plateaus" cen- 
tered at even (a s — ro) and odd (a s = ±oo) values 
of C, respectively. In this regime the underlying step- 
plateau character of the M(a s /ro) function is preserved 
in the M(£) dependence due to the appreciable slope of 
the a s (C) dependence at a s = r and relatively weak di- 
vergence at a s = ±oo. To recapitulate, both for 
and xo 1 and not very deep well the relative quantum 
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defect A^(C,a;o) is an increasing staircase function of the 
depth parameter £ with the steps located at odd {xq <C 1) 
or even (xq » 1) values of £■ The crossover between the 
two regimes can be seen in Fig. 5 as a relatively narrow 
stripe of very weak modulations. 

A qualitatively different staircase-like dependence 
emerges for so > 1 and sufficiently deep well because as 
C — > oo the slope of the a s (() function at a s = r tends 
to zero while the binding resonances, a s — ±oo, become 
progressively more singular. As a result the "steps" and 
"plateaus" switch places - the former becomes centered 
at odd while the latter at even values of (. This is some- 
what similar to what happens in the xq <C 1 limit. This 
observation explains why the small modulation crossover 
stripe in Fig. 5 runs at an angle to the (A^i, () plane. 
The important qualitative difference between the io < 1 
and Xq, C ^> 1 staircases is that the latter have "plateaus" 
centered at half-integer values of M, thus corresponding 
to the coincidence of a node of the function x m Eq.Q 
with the boundary of the inner region r — ro . 



E. Semiclassical treatment 

For most realistic models of the inner potential U s (r) 
the exact analytical calculation of the scattering length 
a s entering the general expression for the quantum defect 
En. lfTTl) may not be possible. Therefore it is pertinent to 
understand whether there is an approximate analytical 
treatment capturing the Zel'dovich spectrum reconstruc- 
tion. This is especially relevant to the xq ^> 1 regime 
when the phenomenon manifests itself only as a modula- 
tion of the quantum defect superimposed on a monotonic 
curve. 

For io > 1 and sufficiently smooth U s (r) the stan- 
dard semiclassical approximation is applicable, and the 
corresponding solution to Eq.(j2j can be written as 



Xsc (r)~(S'(r)r 1/2 sin* 



(15) 



the scattering length, are applicable when the number of 
de Broglie's half- waves So/irh fitting inside U s (r) is very 
large. If we additionally assume the continuity of the 
central potential U(r) in Eq.Q at the boundary of the 



inner region, 



U s (r ) = -h 2 /ma B r , 



(18) 



then En. lTL7Tl simplifies to a s /ro = 1 — (2/xo) ta,n(So/h). 
Combining this with Eq. 114t we find an explicit semiclas- 
sical expression for the quantum defect 



3 
4 



So 

n 



So 

Tit) 



(19) 



where 



which can be interpreted as approximately the sum of 
3/4 — xq/tt, the quantum defect for U s {r) = 0, and the 
number of de Broglie's half- waves So / irh fitting inside the 
inner part of the potential; the estimate fi sc ~ So/irh has 
been given earlier jjj. The number of de Broglie's half- 
waves can be estimated in terms of the dimensionless 
depth of the inner well w ~ mrQ\U s \/h 2 as So/irh~ w 1 ! 2 
which implies that for fixed xq the quantum defect l|19JI 
is a monotonically increasing function of w without any 
modulations. We conclude that the Zel'dovich modula- 
tions of the quantum defect are lost in the semiclassical 
approximation despite the fact that the corresponding 
scattering length itTTjl does exhibit binding resonances. 
Thus for io > 1 a treatment better than semiclassical 
is required to capture the deviations from monotonic be- 
havior; a similar conclusion has been reached earlier [l2T |. 

For the rectangular well of radius r and depth 
Uo characterized by the coupling constant C = 
(SmC/org/V 2 ?! 2 ) 1 / 2 the expression for the scattering 
length i!17|) is exact, and then Ea. (|19fl predicts that 
fj, = 1/4 - xo/tt + C/2 = 1/4 - C/2 < 0. This is the 
xq, C ^> 1 value of the quantum defect in the middle of the 
small modulation stripe in Fig. 5 whose locus, xo = t^C, 
can be deduced from Eci JlSJI . The quantum defect is 
negative because for continuous U(r) the short-distance 
rectangular well potential is always less attractive than 
the Coulomb potential. 



S(r) = / (-2mU s (r)) 1/2 dr 
o 



(16) III. QUANTUM DEFECT OF RYDBERG 

ELECTRON 



is the classical action acquired by a zero-energy particle 
moving radially out from zero to r, and the prime in 
Eq.JEJ denotes differentiation with respect to r. The 
semiclassical expression for the scattering length which 
can be deduced from En ifTBTl with the help of Eq. iflOjl 
has been given by Berry 22] 



Os _ 

ro 



2mr^U s (r ) 



1/2 



. So 
tan — 
h 



(17) 



where S = S(r ). Ens. ltl5)l and ifTTjl . generalizing the 
"rectangular well" expressions for the wave function and 



Now when we understand the manifestations of the 
Zel'dovich effect, and what kind of accuracy is required to 
approximately capture it, we begin computing the quan- 
tum defect of the Rydberg electron as a function of po- 
sition along the Periodic Table. The quantum defect is 
given by the exact result Ea. lfTTll with ro and a s , being 
the size and the scattering length of the residual atomic 
ion, respectively, both dependent upon atomic number 
Z. The resulting dependence will exhibit mod- 

ulations both due to systematic (Zel'dovich) and shell 
effects. As discussed in the Introduction, the shell ef- 
fects obscure systematic trends making it difficult to see 
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that some modulations of (i{Z) have their origin in the 
binding properties of the ionic core. To circumvent this 
inconvenience below we conduct a calculation capturing 
only systematic effects. The comparison of the results 
with both experimental and numerical data (addition- 
ally containing the shell effects) will allow us to disentan- 
gle physically different sources of deviation from purely 
monotonic behavior. 



A. Method of comparison equations 

The short-distance potential U s (r) characterizing the 
residual atomic ion will be assumed to match at its 
boundary the Coulomb potential of unit charge [l0],[12] 
(see Eq.(EJ). As the Rydberg electron moves inside the 
ionic core, the screening of the nuclear charge by the 
inner shell electrons diminishes which implies that for 
r < tq the short-distance potential U s (r) is more attrac- 
tive than the Coulomb potential of unit charge. There- 
fore the quantum defect is a necessarily positive and in- 
creasing function of atomic number Z. As r — ► 0, the 
inner potential approaches that of a nucleus of charge 



Ze, i. e. U s (r — 
Eq.JHl reduces to 



0) - 
fx 

dr 2 



-Ze jr = —Zh/ma B r, and 



2Z 
ra B 



-x = o 



(20) 



This presents a convenient starting point for obtaining 
an approximate solution to the differential equation 
via the method of comparison equations as described by 
Berry and Mount j23|. Since Eq.JSOJ) is exactly solvable, 
and the potentials of Eqs.© and l(20j) are somewhat sim- 
ilar, the solution to © should be also similar to that of 
(1201 and can be transformed into it by a slight deforma- 
tion of coordinates and an amplitude adjustment. The 
details of finding an appropriate mapping are given in 
Ref. [23]; the resulting approximate solution of © is 
then given by 



X(r) 



S(r) 
S'(r) 



1/2 



Jl 



S(r) 



(21) 



The method of comparison equations includes the con- 
ventional semiclassical treatment as a special case [2flj . 
From this more general viewpoint Eq. ITfil can be viewed 
as a result of the deformation and amplitude adjustment 
of the "rectangular well" sine solution. 

To assess the accuracy of (12 H let us first look at the 
limit S(r)/h <C 1. According to Eo. lfTfij) this corresponds 
to r — > when U s (r) — > —Zh 2 /ma B r. Then S(r — > 
0)/H -► (8Zr/a B ) 1/2 and x(r -> 0) ~ r 1 / 2 J^ZZr / a B ) 
which can be recognized as the solution to l(20jl . 

In the opposite limit S{r)/h ^> 1 a semiclassical ap- 
proximation is expected to be valid and Eq. 12H simplifies 

to x(r) ~ (S'(r))~ 1/2 sm(S(r)/h- tt/4). This is similar 
to the naive semiclassical result 1)1511 with the extra phase 



of —tt/4 correcting for the failure of the standard semi- 
classical treatment in the Coulomb field of charge Ze at 
distances r < a B /Z. Thus for the inner potential U s (r) 
which has a Coulombic singularity as r — > but is oth- 
erwise smooth the analog of Ea. lll9Jl is 



(IsC — 



7T 



irk 



(22) 



The expression for the scattering length corresponding 
to Eq.lEU can be found with the help of Eo. lfTo)) 



= 1-4 



ro 




roUjjro) 
U s (r ) 



(23) 

where we also used the condition of continuity itTsl) . 
Eq. jSBt can be used to go beyond the semiclassical ex- 
pression 12211 . Combining Eas. <ll4H and 12. ill we find that 
in the xq,Sq/H ^> 1 limit the quantum defect can be 
presented as fi = fi 8C + Sfi where fi sc is the semiclassical 
answer and the correction, 



Sfi = 



roU' s (r ) 
4nx a U s (r ) 



^ ) - m 



captures the Zel'dovich effect now manifesting itself as a 
simple harmonic modulation superimposed on the semi- 
classical background. The period of the oscillation is ex- 
actly one de Broglie's half-wave while the amplitude is of 
the order Xq . The fact that the latter is independent of 
the number of de Broglie's half- waves fitting inside U s (r) 
implies that the Zel'dovich effect persists for any value of 
Sq/ttH. 



B. Thomas-Fermi model of atomic ion: 
semiclassical solution 

Below we follow Latter ^3] and assume that the po- 
tential of the atomic ion U s (r) can be approximated by 
the Thomas- Fermi theory |24| : 



Us(r) 



h 2 z 
< 

ma B r 



rZ 1 / 3 
ba B 



(25) 



where b = (37r/4) 1 / 2 /2 ~ 0.885, and the universal func- 
tion 4>{y) is the solution to the nonlinear Thomas-Fermi 
equation 



l/2 rf2 

dy 2 



/,3/2 



(26) 



subject to the boundary conditions 0(0) = 1 and </>(oo) = 
[2J|. Then the size of the ion ro and thus the range 
parameter xq are determined by the continuity con- 
dition l|T8jl . i. e. when the Thomas- Fermi potential l(25|l 
meets the Coulomb potential of unit charge: 



,x 2 Z^\ 



1 

z 



(27) 
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Eqs. 12511 and $271 imply that the natural variable 
to characterize the strength of the potential of the 
atomic ion is Z 1 / 3 . Indeed, the typical length scale 
of the Thomas- Fermi theory is ab /Z 1 / 3 , the magni- 
tude of the typical potential is (Zh 2 /raaB){Z 1 / 3 /as) = 
(h 2 /ma 2 B )Z 4 / 3 , and thus the dimensionless coupling con- 
stant w ~ mr1\U s \/h 2 which entered the general analysis 
of Section II is of the order Z 2 I 3 \ the parameter Z 1 / 3 then 
parallels C — w 1 / 2 used in the rectangular well example 
of the inner potential. 

Since the Thomas-Fermi function <fi(y) is a monoton- 
ically decreasing function of its argument, the xq{Z x ^ 3 ) 
dependence defined through Ea J27Jl is a monotonically 
increasing function of Z 1 / 3 . The boundary condition 
0(0) = 1 implies that xq(1) = which is in accor- 
dance with the expectation that for Hydrogen (Z = 1) we 
have the standard Coulomb problem in the whole space 
(xo = 0). As evident from ll27l) small values of <f> are 
relevant for large Z; in view of 4>(y — > oo) — * 144/j/ 3 j24[ 
this means that there is an upper bound to the range 
parameter, x (oo) = 2 13 / 6 3 1 / 3 6 1 / 2 ~ 6.092. 

Another quantity of interest is the number of de 
Broglie's half-waves fitting inside the Thomas-Fermi 
atomic ion 



x 2 Z 1/3 /Sb 



which is a monotonically increasing function of Z 1 / 3 : for 
Z — ► 1 it vanishes as xq/it while for large Z we have 
Sq/ttH ~ Z 1 ' 3 , a known result 0. 

For intermediate values of the atomic number the Z 1 / 3 
dependences of the range parameter xq and the number 
of de Broglie's half- waves So/irh can be found by numer- 
ically solving the Thomas- Fermi equation 12 (ill , inverting 
(l27l) . and computing the integral ffify . The results are 
displayed in Fig. 6 where we show xo/n and So/nh as 
functions of Z 1 / 3 . These functions are used to also plot 
the semiclassical quantum defect given by Eq. 122t . To as- 
sess the accuracy of the resulting /i sc (Z 1 / 3 ) dependence 
we need to verify whether approximations used to derive 
Ea, l|22H are adequate. The first assumption, Xq ^> 1, is 
equivalent to the assertion that the quantum defect for 
U s (r) — can be replaced by its large x$ limit. Looking 
at Fig. 1 where the dependences in question are compared 
we conclude that "large" here really means xq > 7T. The 
second, semiclassical assumption Sq/H 3> 1 in practice 
has a good accuracy provided So /irh, the number of de 
Broglie's half- waves fitting into the inner potential U 8 {r), 
is anything more than one or two. Inspecting Fig. 6 we 
see that the conditions xq/-k > 1 and So/irK > 2 are 
satisfied for Z > 8. This is also the practical condition 
for the Thomas-Fermi-Latter model of the residual ion 
to be applicable. We note additionally that although 
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z 1/3 

Figure 6: Plots of the number of de Broglie's half-waves 
So/nh, 128II . the range parameter xo/n, 127II . and semiclas- 
sical quantum defect fi S c, i'2'2l as functions of Z 1 ^ 3 for the 
Thomas- Fermi model of the residual atomic ion. 

Eq.Q2l is not expected to be valid for smallest Z, the 
limit n(l) — is nevertheless correctly reproduced and 
for any Z the semiclassical quantum defect is an increas- 
ing positive function of Z 1 ! 3 in accordance with physical 
expectation. We conclude that except possibly for the 
elements of the first row of the Periodic Table, the semi- 
classical result for the quantum defect En. 122^ shown in 
Fig. 6 is accurate; only qualitative agreement is expected 
for lightest elements. 

C. Beyond semiclassical approximation: connection 
to binding properties of ionic core 

A more accurate ^(Z 1 / 3 ) dependence can be found 
by computing the scattering length ll23ll and substitut- 
ing the outcome together with the xq{Z x ^ 3 ) dependence, 
Fig. 6, in our general expression for the quantum defect 
lUTf . The result is shown in Fig. 7 where we also plot the 
semiclassical quantum defect, ^ sc , Eq.J22J (gray scale). 
It now becomes obvious that for any Z the bulk contribu- 
tion into the quantum defect is well-captured semiclassi- 
cally. The Zel'dovich spectral modulation clearly visible 
in Fig. 7 is a relatively weak effect. To separate the mod- 
ulation from the monotonic s semiclassical background 
the inset shows the difference Sfi = /t — jjL sc which ap- 
pears to be a nearly periodic function of Z 1 / 3 . To better 
understand the meaning of this periodicity the inset also 
shows the limiting expression 12411 (gray scale) 25]. Both 
curves are exactly in phase and for sufficiently heavy el- 
ements their magnitudes agree semi-quantitatively. This 
observation implies that the Zel'dovich modulation is a 
periodic function of the number of de Broglie's half- waves 
fitting inside the ionic core of the atom. 

A complementary way to see the connection between 
the spectral modulation and the binding properties is 
presented in Fig. 8 where we compare the <i\x{Z x l 3 ) de- 
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z 1/3 

Figure 7: Dependence of the quantum defect /i on Z X ^ A along 
with its semiclassical approximant /i sc , Ea. l22ll (gray scale). 
The inset shows the Zel'dovich modulation S/i(Z 1 /' 1 ) = /i — /i sc 
together with the limiting expression, En. 12411 (gray scale). 
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Figure 8: Dependences of the Zel'dovich modulation S/j, = 
fi — /Use and the reduced scattering length of the ionic core 
Os/ro (gray scale) on Z 1 ^ 3 . The lines a s /ro — and a 3 /ro = 1 
are also shown to help the eye. 



pendence with the behavior of the scattering length of 
the residual atomic ion II23H . The latter, numerically 
computed for the Thomas-Fermi model of the residual 
atomic ion, Eqs.JlSJ - is shown in gray scale. The 

binding singularities of the scattering length are nearly 
equidistant confirming the earlier observation that the 
parameter Z 1 / 3 is analogous to ( used in Figs. 3 and 5 
to display the Zel'dovich effect for the rectangular well 
model of the inner potential. A qualitatively similar be- 
havior of the scattering length of the Thomas- Fermi atom 
as a function of Z has been reported by Robinson [2flj : 
quantitative differences may be attributed to the assump- 
tion j2fj that the Thomas-Fermi potential vanishes at a 
distance of the order clb/Z 1 ^ which is different from our 
choice of r , Eas. lfLsj) and (23 . 

Fig. 8 makes it clear that the maxima of the oscillation 
5/j, occur when a s — rn, i. e. when the node of the 
function x( r ) m Eq.lQJ coincides with the boundary of 
the atomic ion. On the other hand, the minima of Sfi 
are correlated with binding singularities of the ionic core, 
a s = ±oo, thus corresponding to the antinode of x(r) 
being near the ion boundary. 



IV. COMPARISON WITH EXPERIMENTAL 
AND NUMERICAL DATA 

We found experimental values of quantum defects for 
37 elements of the Periodic Table. These data and their 
sources are compiled in Table 1 where we also list sys- 
tematic quantum defects of our work (also displayed in 
Fig. 7) . Some of the figures which we re gard as "exper- 
imental" came from an on-line database 29] where the 
quantum defect is computed from available spectroscopic 
data. 

In cases of He, Be, Mg, Ca, and Mo there is more than 
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Figure 9: Systematic, experimental, and numerical depen- 
dences of the quantum defect /i on Z 1 ' 3 . To help orientation 
within the Periodic Table experimental alkali data are circled. 

one value of the quantum defect available depending on 
the angular momentum of the ionic core of the atom. 
Since our theory represents average properties and does 
not distinguish between different LS terms of an atomic 
configuration, in Table 1 we chose to show only the values 
corresponding to lowest angular momentum of the ionic 
core. It turns out they better agree with our calculation 
than those left out. 

The works of Manson and Fano, Theodosiou and 
Dehmer 0] contain graphs of numerically evaluated 
)jl{Z) dependence for all elements. After verifying that 
the results of both studies are nearly identical, we chose 
to restrict ourselves to those of the later Ref.frj]. 

Experimental, numerical and systematic ^(Z 1 / 3 ) de- 
pendences are displayed in Fig. 9. In order to produce 
the numerical curve, the data 0] have been scanned, 
digitized, and replotted as a function of Z 1 / 3 . We also 
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Table I: Experimentally measured quantum defects for series 
of elements with their atomic numbers Z and corresponding 
references. Systematic quantum defects of this work are also 
displayed for comparison. 



z 


Element 


Experimental fi 


Reference 


Systematic fi 


2 


He 


.139 


127] 


.110 


3 


Li 


.400 


Tool 

I28| 
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circled locations of alkali metals because their ionic cores 
have noble element electronic configurations thus mark- 
ing (for the Rydberg atom problem) the end of a period. 

Fig. 9 makes it clear that all three dependences are 
in fairly good agreement and Z 1 / 3 is certainly the right 
variable to use for analysis. It is not surprising that nu- 
merical results ^3] are generally in better agreement with 
experimental data than our systematic findings because 
our calculation omits the shell effects. 



A. Effects of shell structure 

In order to be able to separate systematic and shell ef- 
fects, in Fig. 10 we display <5/i = /i — fi sc , the modulation 
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Figure 10: Modulation of systematic, experimental and nu- 
merical quantum defect <5/i = /i — /i sc relative to the semiclas- 
sical background, En. 12211 as a function of Z 1 ^ 3 . 



of the quantum defect relative to the monotonic semiclas- 
sical background, Eci. il22|) . which accounts for the bulk of 
the quantum defect value. Fig. 10 shows that the exper- 
imental and numerical variations of the quantum defect 
are bounded which is consistent with the view that they 
are due to repetitive physics. Moreover, the systematic 
modulation due to the Zel'dovich effect appears to have 
an amplitude which is several times smaller than those 
of experimental and numerical data. This observation 
implies that it may be possible to understand gross fea- 
tures of the experimental and numerical modulations of 
the quantum defect as mostly due to the effects of the 
shell structure. 

This viewpoint can be supported by qualitative anal- 
ysis which rests on the semiclassical result 12211 . First, 
let us anticipate the outcome of incorporating the shell 
effects into the calculation. This amounts to replacing 
the smooth inner potential U s {r) by one with modula- 
tions due to the spatial variation of the electron density 
reflecting the shell structure. This replacement will re- 
sult in a value of the size of the ionic core r , Ea. fTsl) . 
generally different from its systematic counterpart. 

Let us additionally assume that the inner potential 
U s (r) with shell effects included is still sufficiently smooth 
so that a semiclassical treatment is valid. The corre- 
sponding quantum defect (l22l) will deviate away from the 
systematic result due to different values of the range pa- 
rameter xq, Eq.Q, and the number of de Broglie's half- 
waves, So/ttH. Because the latter involves the integral of 
(-[/s(r)) 1 / 2 from zero to r (see Eqs.HEl and <GBJl) the 
modulations above and below the systematics present in 
U s (r) are expected to largely cancel each other and the 
deviation from our results can be mostly attributed to 
the different size of the ionic core. 

This argument implies that the quantum defect is 
strongly sensitive to the value of the size of ionic core 
ro and weakly sensitive to the details of the inner po- 
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Figure 11: Slater's ionic radii for singly-charged positive ions 
together with a series of corresponding ionic radii in crystals 
and systematic sizes of ionic core of the Rydberg atom, all in 
atomic units, as functions of Z 1 ^ 3 . Numerical variation of the 
quantum defect 5fX = fJ, — [i S c is also displayed to show the 
correlation with Slater's radii. 

tential U s (r). In reality the inner potential may not be 
smooth enough for the semiclassial treatment to be quan- 
titatively correct. Therefore we do not expect more than 
a qualitative insight into the trends of the variations of 
the quantum defect induced by the shell effects. 

The simple rule that emerges can be most easily de- 
duced from Fig. 6 by keeping in mind the relationship 
between the range parameter xq, Eq.Q and the size of 
the ionic core ro oc Xq : deviation in tq away from system- 
atics leads to the same sign deviation in the quantum de- 
fect. Since the size of the ionic core of the Rydberg atom 
has a physical meaning close to that of an ionic radius, 
to verify the correlation we need a set of ionic radii for 
singly-charged positive ions as a function of position Z 
along the Periodic Table. 

Seventy five years ago J. C. Slater gave a very use- 
ful, general, empirical set of rules to approximate analyt- 
ically atomic wave functions for all the elements in any 
stage of ionization. The radial part of the single-electron 
wave function is selected in the form 

tP(r) ozr n *- 1 e- z * r / n ' aB (29) 

which can be recognized as the large-distance asympotics 
of a Hydrogen-like wave function with an effective quan- 
tum number n* and an effective nuclear charge Z*e. 
Based on the underlying electronic structure, Slater's 
rules assign values of n* and Z* to the electrons of each 
shell of an atom or ion, so that a complete set of single- 
electron wave functions can be constructed. 

For a given shell the maximum of the electron density 
4:7rr 2 ip 2 (r) is located at 

r ma x/a B = {n*) 2 /Z* (30) 

which formally coincides with the expression for the ra- 
dius of the corresponding circular orbit in Bohr's old the- 



ory. The radius of the maximum density of the outermost 
shell is expected to correlate with the size of the atom or 
ion. Specifically, Slater defines an ionic radius r > r max 
as a distance at which the electron density becomes 10% 
of its maximal value. 

Using the existing knowledge of electronic configura- 
tions j43| we applied Slater's rules to calculate the ionic 
radii of singly-charged positive ions. The result is shown 
in Fig. 11 where the elements marking the beginning or 
an end of more dramatic changes in the ionic radius are 
labeled. For comparison we also displayed a series of 
ionic crystal radii [4J| used for predicting and visualizing 
crystal structures. Crystal ionic radii are based on exper- 
imental crystal structure determinations, empirical rela- 
tionships, and theoretical calculations. As Fig. 11 shows, 
they are in fair agreement with their Slater's counter- 
parts. We hasten to mention that neither Slater's nor 
the crystal ionic radii are expected to coincide with what 
we define as the size of the ionic core of the Rydberg 
atom, Eq. flU - It seems highly plausible, however, that 
Slater's ionic radii are correlated with the sizes of ionic 
core of the Rydberg atom. 

Inspection of Fig. 11 tells us that the average Slater's 
ionic radius slowly grows with Z 1 / 3 in fairly good agree- 
ment with our systematic result. A closer look reveals 
that our systematic radius appears to be consistently 
smaller than its Slater's counterpart. If the same re- 
lationship would hold between the systematic and (un- 
known) exact sizes of ionic cores of Rydberg atoms, then 
the fact that experimental and numerical quantum de- 
fects in Figs. 9 and 10 are generally larger than their 
systematic counterparts would be explained. 

The large variation of the Slater's ionic radius away 
from the average trend is due to the effects of the shell 
structure. Their role in determining the ionic radius can 
be most easily visualized based on the expression for the 
radius of the maximum electron density 13011 which cor- 
relates with the ionic radius. This result emphasizes the 
following main principles: 

(i) As Z increases, all the n levels move down in energy 
which amounts to replacing n by its effective counterpart 
n* < n. If the electrons are added to an outer shell, the 
effective nuclear charge Z*e seen by each of them grad- 
ually increases. This is because the outer shell electrons 
are relatively inefficient in shielding the nuclear charge. 
As a result the ion slowly contracts. 

(ii) As a new outer shell begins to fill, the effect of going 
into the higher shell outweighs the effect of lowering of 
an n level as Z increases to Z + l. This corresponds to an 
abrupt increase of the effective principal quantum num- 
ber n*. Moreover, the effective charge Z*e seen by the 
outermost electron drops because now all remaining Z— 2 
electrons belong to inner shells thus efficiently screening 
the nuclear charge. These changes in n* and Z* cause a 
sharp increase of the ion size. 

Slater's rules [42] add a quantitative aspect to these 
principles. In the following explanation of the variation 
of the ionic radius, Fig. 11, we are always speaking of the 
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positive singly-charged ions whose electronic configura- 
tions are taken from the NIST database 

In going from 2 He (Is 1 ), to 3 Li (Is 2 ) the ion size 
decreases. As one moves to 4 Be ([Heps 1 ) the added 
electron enters the higher 2s shell, and the ionic ra- 
dius sharply increases. Similar increases take place 
as one goes from every alkali to the following alkali 
earth ion. A related jump in ionic radius also occurs 
past every noble element ion, 29 Cu ([Ar]3d 10 ) — > 30 Zn 
([Ar^Hs 1 ), 47 Ag ([Kr]4d 10 ) -> 48 Cd ([Kr]4d 10 5s 1 ), 
and 79 Au ([Xe]4/ 14 5d 10 ) -» 80 Hg ([Xe]4/ 14 5d 10 6s 1 ) be- 
cause a higher s-shell starts to be occupied. An analogous 
argument explains the sharp increase of ionic radius while 
going from 24 Cr ([Ar]3d 5 ) to 25 Mn ([Ar^dMs 1 ) and from 
42 Mo ([Kr]4d 5 ) to 43 Tc ([Kr]4d 5 5s 1 ) 

As one moves from 4 Be ([He]2s 1 ) to n Na ([He]2s 2 2p 6 ) 
the 2s and 2p shells are filled by the electrons, and the 
ion size gradually decreases. A similar effect explains the 
decrease of the ionic radius along the 12 Mg ([Ne]3s 1 )- 19 K 
([Ne]3s 2 3p 6 ), 30 Zn ([Ar]3d 10 4s 1 )- 37 Rb ([Ar]3d 10 4s 2 4p 6 ), 
48 Cd ([Kr]4d 10 5« 1 ) JS6 Cs ([Kr]4d 10 5s 2 5p 6 ), and 80 Hg 
([Xe]4/ 14 5d 10 6s 1 )- 87 Fr ([Xe]4/ 14 5d 10 6s 2 6p 6 ) sequences. 

The ionic radius decreases through the 20 Ca ([Ar]4s 1 )- 
21 Sc ([Ar]3d 1 4s 1 )~ 22 Ti ([Ar]3d 2 4s 1 ) segment. This hap- 
pens because the electrons filling the 3d shell only par- 
tially screen the nuclear charge - as a result the outer 
4s electron sees a gradual increase of effective Z* . The 
same argument explains the decrease of ion size while go- 
ing from 25 Mn ([Ar]3d 5 4s 4 ) to 26 Fe ([Ar^d^s 1 ) which 
is merely a continuation of the Ca-Ti segment. The de- 
crease of ionic radius through the 73 Ta ([Xe]4/ 14 5d 3 6s 1 )- 
77 Ir ([Xe]4/ 14 5c? 7 6s 1 ) series can be similarly under- 
stood. In fact, the first entry of this series is 70 Yb 
([Xe]4/ 14 5d°6s 1 ) where we intentionally modified the 
standard notation to show the absence of the 5d elec- 
tron. 

The size of the ion first abruptly decreases while 
going from 22 Ti ([Ar]3d 2 4s 1 ) to 23 V ([Ar]3d 4 ) and 
then continues decreasing more gradually as one moves 
to 24 Cr ([Ar]3d 5 ). The sudden change is due to the 
fact that the outer shell changes from 4s to 3d which 
can be viewed as a decrease in the effective quantum 
number n* . The subsequent slower increase of the ionic 
radius is due to the increase of the effective nuclear 
charge Z*e seen by the larger number of 3d electrons. 
The same trend is exhibited in the 40 Zr ([Kr]4d 2 5s 4 )- 
41 Nb ([Kr]4d 4 )- 42 Mo ([Kr]4d 5 ) sequence. A very 
similar behavior is found in the 26 Fe ([Ar]3d 6 4s 1 )- 
27 Co ([Ar]3d 8 )- 28 Ni ([Ar]3d 9 )- 29 Cu ([Ar]3d 10 ), 43 Tc 
([Kr]4d 5 5s 1 )- 44 Ru ([Kr]4d 7 )- 45 Rh ([Kr]4d 8 )- 46 Pd 
([Kr]4d 9 )- 47 Ag ([Kr]4d 10 ), and 77 Ir ([Xe]4/ 14 5d 7 6s 1 )- 
78 Pt ([Xe]4/ 14 5d 9 )- 79 Au ([Xe]4/ 14 5d 10 ) series. 

Superficially, a similar steep decrease of the 
ion size is followed by more gradual decrease in 
the 70 Yb ([Xe]4/ 14 6s 1 )- 71 Lu ([Xe]4/ 14 6s 2 )- 72 Hf 
([Xe]4/ 14 5d 4 6s 2 ) sequence. This behavior can be ex- 
plained by noticing that in both steps the effective Z* 
felt by a 6s electron increases; the increase during the 



second step is smaller because inner shell d electrons are 
more efficient in screening the nuclear charge than s elec- 
trons. 

While going from 38 Sr ([Kr]5s 4 ) to 39 Y ([Kr]5s 2 ) and 
to 40 Zr ([Kr]4d 2 5s 1 ) a decrease of ionic radius follows by 
an increase. This happens because the effective Z* felt 
by an outer 5s electron first increases and then decreases. 
The increase of the size of the ion is somewhat smaller 
than the decrease because two Ad electrons in Zr only 
partially shield two extra units of the nuclear charge. 

For the 63 Eu ([Xe]4/ 7 6s 1 )- 64 Gd ([Xe]4/ 7 5d 1 6s 1 )- 
65 Tb ([Xe]4/ 6s ) sequence the ionic radius first de- 
creases and then increases back to its initial value. This 
happens because the d electron in Gd is less effective in 
shielding the nuclear charge than the / electron in Tb 
(considered perfect in Slater's scheme). 

One of the less intuitive increases of ionic radius takes 
place while going from 56 Ba ([Xe^s 1 ) to 57 La ([Xe]5d 2 ). 
On one hand, the effective principal quantum number 
n* decreases which according to il3f)|l should lower the 
ionic radius. However as compared to the 6s electron 
of Ba, the effective Z* seen by one of La's 5d electrons 
also decreases. This happens because the inner shell 5sp 
electrons screen the nuclear charge more effectively if the 
outer shell electrons are in a d state (La) as compared to 
an s state (Ba). As a result the decrease in Z* outweighs 
the decrease in n* thus leading to an increase of the ion 
size. A very similar argument explains the reversed de- 
crease of ionic radius taking place as one goes from 58 Ce 
([Xe]4/ 1 5d 2 ) (whose size is identical to that of La) to 
59 Pr ([Xe]4/ 3 6s 1 ) (identical in size to Ba). Here the role 
of the 4/ electrons merely reduces to compensating for 
the increase of the nuclear charge. 

The ionic radius does not change as one moves from 
57 La ([Xe]5d 2 ) to 58 Ce ([Xe]4/ 1 5d 2 ) because the increase 
of nuclear charge is exactly compensated by adding a 
4/ electron. The same argument explains the con- 
stancy of the ionic radius along most of the lanthanide 
sequence, 59 Pr ([Xe]4/ 3 6s 1 )- 63 Eu ([Xe]4/ 7 6s 1 ), 65 Tb 
([Xe]4/ 9 6s 1 )- 70 Yb ([Xe]4/ 14 6s 1 ). The only exception 
from this trend, 64 Gd ([Xe]4/'5d 6s 1 , has already been 
discussed. 

Now when the variations of the Slater's ionic radius 
are understood, we can compare them with experimental 
and numerical modulations of the quantum defect. Since 
the numerical data are more extensive than experimen- 
tal findings, and the agreement between the two is fairly 
good, in Fig. 11 we only show the variation of the nu- 
merically evaluated quantum defect 5/j, = fj, — jj, sc . The 
inspection of Fig. 11 leaves no doubt that the variations 
of the quantum defect with Z 1 / 3 are correlated with those 
of the Slater's ionic radius - even the minute changes of 
the latter find their way in the corresponding changes of 
the former. 

There are however two places where it appears there 
is a disagreement with expectation: 

(i) Along most of the lanthanide sequence the Slater's 
ionic radius does not change while the variation of the 
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quantum defect decreases. This can be understood as an 
artifact of the Slater's rules. In reality the / electrons 
do not perfectly screen the nuclear charge. Therefore the 
effective charge seen by the outer 6s electron increases 
with Z 1 / 3 and correspondingly the ionic radius should 
decrease. Then the quantum defect variation should de- 
crease as well which is in correspondence with numerical 
results. We also note that the Gd dip of the ionic radius 
is reproduced in the quantum defect variation. 

(ii) For the elements past Hg, the Slater's ionic radius 
decreases with Z 1 / 3 while the variation of the quantum 
defect increases. The reason why it happens is unclear. 
It cannot be ruled out that here the effects of the shell 
structure might be so strong that our correlation rule 
derived from semiclassical arguments breaks down qual- 
itatively. 

Overall, the analysis of this section makes it certain 
that the gross features of the quantum defect variation 
with Z 1 / 3 are due to the effects of the shell structure. 



B. Zel'dovich modulation 

Qualitative analysis is of little use in trying to see the 
Zel'dovich effect in experimental and numerical data be- 
cause the Zel'dovich modulation has an amplitude which 
is several times smaller than that due to the effects of the 
shell structure (see Fig. 10). A way to proceed quantita- 
tively is suggested by the fashion in which the shell and 
systematic effects are coupled. 

Inspecting the limiting expression for the Zel'dovich 
modulation, Eci. il24|) . it is straightforward to see that af- 
ter replacing systematic U s (r) with the one accounting 
for the shell structure, the amplitude of the Zel'dovich 
modulation will become strongly sensitive to the effects 
of the shell structure because it is determined by the log- 
arithmic derivative of the inner potential U s (r) and the 
range parameter xq. On the other hand, the period of 
the oscillation is far less sensitive to the effects of the 
shell structure since it is determined by the number of 
de Broglie's half-waves fitting inside the ionic core of the 
atom. 

This last observation suggests that it might be possible 
to see the Zel'dovich effect in the experimental and nu- 
merical Fourier spectra of the quantum defect variation 
5(1 = (i— fx sc (see Fig. 10) as a peak whose location can 
be brought in correspondence with the systematic theory. 
To proceed in this direction, in a range of Z 1 / 3 of length 
L we expand the quantum defect variation into a Fourier 
series 

S^Z 1 ' 3 )^ Mfc cxp^/cZ 1 / 3 ) (31) 

k=2vp/L 

where p = 0,±l,±2,.... In order to numerically evaluate 
the Fourier coefficients /ife = the experimental and 
numerical S^Z 1 / 3 ) dependences (see Fig. 10) were fitted 
with a cubic spline which was then sampled equidistantly 
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Figure 12: Systematic, experimental and numerical ampli- 
tudes of the Fourier coefficients of the quantum defect varia- 
tion (arbitrary units, same normalization ) \jjbu \ as functions of 
k for k > 0. The peaks at k ~ 11 correspond to the Zel'dovich 
effect. 

in Z 1 ! 3 to extract the Fourier spectrum. The result for 
the magnitude of the Fourier coefficients \\ik\ as a func- 
tion of k is displayed in Fig. 12 as a series of solid dots 
which for convenience are connected by straight line seg- 
ments. The uncertainty of the location of each dot along 
the k axis, 2ir/L, is the distance between the nearest 
values of k. Since the last available experimental quan- 
tum defect corresponds to 83 Bi (see Table 1), for both 
experimental and numerical data we restricted ourselves 
to L = 83 1 / 3 . The finite values of fio correspond to the 
presence of nonzero background in experimental and nu- 
merical Sfi^Z 1 / 3 ) dependences and are of no interest to 
us. 

For comparison in Fig. 12 we also show the Fourier 
spectrum of our systematic calculation which, as ex- 
pected, has only one peak corresponding to the Zel'dovich 
effect. The position of the peak along the k axis can be 
understood from the large Z asymptotics of the Thomas- 
Fermi action Sq/H ~ 5.2Z 1 / 3 l|28jl. Comparing this with 
the limiting expression for the Zel'dovich modulation 12411 
we would expect a peak at k ~ 10.4. The peak in 
Fig. 12 is located at a slightly different value of k ~ 11 
which is due to the fact that the asymptotic behavior 
So/H ~ 5.2Z 1 / 3 becomes numerically accurate only for 
Z 1 / 3 exceeding 10. 

Both the experimental and numerical spectra in Fig. 12 
have peaks at the same value of k ~ 11 which we argue 
are the signatures of the Zel'dovich effect. 

It is curious that both experimental and numerical 
spectra have another peak in common located at k ~ 5.5. 
This peak which is about twice as high as that due to the 
Zel'dovich effect is natural to relate to the effects of the 
shell structure. The existence of this peak translates into 
the (2tt/5.5)Z 1 / 3 = 1.14Z 1 / 3 periodicity of the quantum 
defect variation due to the effects of the shell structure. 
This conclusion resembles the Z 1 ^ 3 periodic oscillation of 
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the ground-state energy of an atom away from the sys- 
tematic trend ^]|. With uncertainty of the peak location 
in mind, one may speculate that our result is a manifes- 
tation of the same effect for highly-excited states. More 
work is necessary to bring understanding to this issue. 

We also repeated the same Fourier analysis by choos- 
ing the range of Z 1 ' 3 to be L = 102 1 / 3 which includes all 
the numerical data. As far as experimental data go in the 
range of Z 1 / 3 between 83 1 / 3 and 102 1//3 we used extrap- 
olation of our cubic spline fit. As a result the peaks just 
discussed slightly change their positions and amplitudes 
but within the 2ir/L uncertainty systematic, experimen- 
tal and numerical Fourier spectra share a peak in com- 
mon corresponding to the Zel'dovich effect. Similarly, 
experimental and numerical spectra continue to share a 
peak in common due to the effects of the shell structure. 



V. CONCLUSIONS AND FUTURE 
DIRECTIONS 

In this paper we analyzed in a model-independent fash- 
ion the weakly-bound s spectra of the distorted Coulomb 
problem for arbitrary relationship between the range of 
the inner potential and Bohr's radius of the Coulomb 
field. We demonstrated that the spectra are fairly sensi- 
tive to the binding properties of the inner potential which 
constitutes the essence of the Zel'dovich effect, and es- 
tablished the corresponding details of spectral changes. 
Armed with these results, we conducted an analysis of 
experimental and numerical Rydberg spectra of atoms 
along the Periodic Table which indeed show an evidence 
of the Zel'dovich effect. Our analysis can be extended 
and adopted in several directions: 

First, there is an abundance of experimental and nu- 
merical data for atomic Rydberg states of finite angular 



momenta which are likely to contain signatures of the 
Zel'dovich effect. However, in the limit of a very short- 
ranged inner potential the way the effect manifests it- 
self is somewhat different from its s state counterpart 
0. This observation makes it pertinent to generalize our 
analysis to the case of finite angular momentum. 

It has been known for some time .4.6] that the Ryd- 
berg formula Q is superior to the Wannier (Bohr) for- 
mula quoted in textbooks Q in representing excitonic 
spectra in condensed matter systems. Experimental ex- 
amples here include clean and doped rare-gas solids and 
rare-gas impurities in solid hydrogen. Although this is 
a context in which the Zel'dovich effect has been origi- 
nally discovered [3], to the best of our knowledge there 
were no attempts to relate it to excitonic quantum de- 
fects. Because of the dielectric screening of the Coulomb 
interaction, the Zel'dovich effect in these systems is ex- 
pected to be more pronounced than in atomic Rydberg 
spectra. Only minor changes to our analysis are needed 
to understand the excitonic Rydberg spectra. 

Other examples of systems where the Zel'dovich ef- 
fect should have experimental signatures include Ryd- 
berg ions and electronic image states 
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